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Abstract. We consider the electronic transport properties of phosphorus (P) doped silicon nanowires (SiNWs). 
By combining ab initio density functional theory (DFT) calculations with a recursive Green's function method, 
we calculate the conductance distribution of up to 200 nm long SiNWs with different distributions of P dopant 
impurities. We find that the radial distribution of the dopants influences the conductance properties significantly: 
Surface doped wires have longer mean-free paths and smaller sample-to-sample fluctuations in the cross-over from 
ballistic to diffusive transport. These findings can be quantitatively predicted in terms of the scattering properties 
of the single dopant atoms, implying that relatively simple calculations are sufficient in practical device modeling 
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1. Introduction 

The continuous reductions in feature sizes in the semi- 
conductor industry have generated much interest in 
novel one-dimensional nanostructures such as silicon 
nanowires (SiNWs) [TJ [2]. Scattering by defects or 
dopants will be increasingly important with decreas- 
ing sizes, and at the same time sample-to-sample vari- 
ations become a crucial issue: when the device length 
and the mean free path are comparable, and shorter 
than the coherence length, variations of the positions 
of the individual dopant atoms can affect the conduc- 
tance of the wire significantly. Accurate models for the 
electronic transport properties of nanowires including 
sample- variations are clearly desirable. 

The mathematical theory of the conductance of dis- 
ordered quasi one-dimensional systems has reached a 
high level of understanding in the diffusive as well as 
the localization regime [3 0] . The diffusive regime is 
characterized by l e < L < £ with L, l e and £ being the 
sample length, mean free path (MFP) and localization 
length, while in the localization regime L > £. The 
universal conductance fluctuations (UCF) is a remark- 
able property of disordered systems: In the diffusive 



regime the sample-to-sample variations are described 
by a Gaussian distribution with universal width of the 
order of e 2 /h independent on the shape of the conduc- 
tor and type of disorder [5] . 

Several recent theoretical works used density- 
functional theory (DFT) to consider the energetics of 
dopants in SiNWs [6j [Tj . Further, Fernandez-Serra et 
al.[S] considered scattering properties of single phos- 
phorus (P) dopants. In a recent work [5] we studied 
sample-averaged conductance properties ((G), std(G), 
MFP) in long SiNWs containing a random distribu- 
tion of either B or P dopants along the wire. In this 
work we continue the analysis focusing more on the 
conductance distribution and on the influence of the 
radial dopant profile. 

2. Method 

In this work we model the conductance properties of 
SiNWs with randomly distributed P dopant atoms. 
The wires are oriented along the [1 0] direction and 
are surface passivatcd by hydrogen atoms to avoid dan- 
gling bonds. The Hamiltonians describing the pristine 
wire as well as regions containing dopant impurities 
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are obtained using ab initio density-functional theory 
(DFT) calculations [TT]. The DFT calculations on the 
doped regions use supercells containing nine wire unit 
cells (837 atoms) with a total length of 50.4 A and 
with the single dopant atom placed in the middlemost 
unit cell. The atomic positions of all the atoms in the 
central unit cell containing the dopant atom as well 
as the two neighboring unit cells have been fully re- 
laxed, until the maximum force was smaller than 0.04 
eV/A. The large super cell ensures that the there is 
no dopant-dopant interaction and that the electronic 
structure at the supercell-supercell interface is very 
close to that of the pristine wire. The length and en- 
ergy dependent conductance of a long wire with ran- 
domly placed P atoms is calculated using a standard 
recursive Green's function (GF) method [TO]. In each 
recursion step, we add a unit cell from cither the pris- 
tine wire or from a region around a dopant atom. At 
each energy we do this calculation for 560 different 
realizations of the dopant positions to obtain the con- 
ductance distribution. By changing the relative occur- 
rence of different dopant positions, we can model dif- 
ferent radial dopant profiles. In all calculations we use 
an average dopant-dopant separation d = 10 nm corre- 
sponding to a bulk doping density of n « 10 19 cm~ 3 . 
Since P dopants arc n-type we consider only energies 
in the conduction band and measures all energies rel- 
ative to the conduction band edge, E c . 

3. Single-dopant results and estimates 

Fig. [1] shows the transmission vs. energy through in- 
finite wires containing only a single P dopant atom 
placed at five different substitutional positions (1-5) 
indicated in the inset. The P dopants located in 
the bulk positions (1-3) generally scatter the electrons 
more than the surface positions (4-5). 

In a recent work [9] we showed that the sample- 
averaged properties, (G), std(G) and MFP could be 
accurately predicted using information only from the 
single-dopant transmissions shown in Fig.[TJ We define 
an average scattering resistance 
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position into account. In the present study i = 1 ... 5, 
corresponding to the five different dopant positions 
shown in the inset in Fig. [1] Denoting the average 
dopant-dopant distance, d, the average resistance of a 
wire of length, L, is estimated as 



R(E, L) = R C (E) + (R S (E)} L/d. 



(2) 




Figure 1: Single dopant conductance vs. energy. Each curve 
correspond to a specific position of the dopant atom, as in- 
dicated in the inset showing a cross section of the wire with 
the possible dopant positions marked. 



From the linear region of the (R) vs. L curves (see 
Fig. [3] (a)) one can extract the energy dependent MFP, 
l e (E) through the relation 



(R(E, L)) = R C (E) + R C {E) L/l e {E). 



(3) 



Comparing with ^ we obtain an expression for the 
MFP given only by the single-dopant transmissions, 
the doping-dopant distance, and the dopant distribu- 
tion: 
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We emphasize that the relations above only hold for 
wire lengths shorter than the localization length, £. 
We can calculate £ using the relation [3], 

Z{E) = \{N{E) + l)l e {E), (5) 

which now can be estimated using the single-dopant 
result, l e (E) [9]. 

4. Long wire results 



where E^ift = 1 and 1/G (E) = h/(N(E)2e 2 ) = 
R C (E) is the contact resistance with N(E) being the 
number of conducting channels at energy E. The first 
term on the right is the average of the single-dopant 
conductances taking the probability, p,-, of each dopant 



Fig. [5] shows the conductance distribution for differ- 
ent wire lengths at energy E = 0.16 cV (a)-(d) and 
E = 0.01 eV (e)-(h). The radial dopant distribution is 
uniform, i.e. the positions 1-5 are equally likely. The 
conductance distribution at the higher energy (left) 
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develops from a quasi-ballistic regime (a) with a few 
peaks in the histogram, to the diffusive regime (b)+(c) 
and finally to the beginning of the localization regime 
(d). 

In the diffusive regime the conductance distribu- 
tion is well described by a Gaussian with the stan- 
dard deviation being close to the universal conduc- 
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tance fluctuation (UCF) value of 0.37 • -=|-. In the 
ballistic and diffusive regimes, the sample averaged 
resistance (R) = I /(G) increases linearly with wire 
length, as illustrated in Fig. [3] (a) (circles). For wire 
length L > 150 nm the resistance starts to increase 
above the initial linear behavior, thus entering the lo- 
calization regime. Notice that the localization sets in 
when (R) rs 1 • in agreement with Thouless' ar- 
gumentation [T2] . The resistance of the surface-doped 
wires Fig. [3] (a) (stars) increases linearly throughout 
the length range and stays in diffusive regime. The two 
solid lines in Fig. [3] (a) obtained using © clearly re- 
semble the sample-averaged values (circles and stars) 
in the linear regimes. This demonstrates, that the av- 
erage resistance in the ballistic- and diffusive regimes 
can be accurately predicted for various dopant dis- 
tributions only based on the single-dopant scattering 
properties. 
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Figure 2: Conductance distributions at different wire lengths 
(rows) for a uniform radial dopant distribution at energy 
E = 0.16 eV (left) and E = 0.01 eV (right). 

In the distribution plots the diffusive-localization tran- 
sition starts when the low conductance tails of the 
histograms come close to zero and starts to 'pile up'. 
Eventually the distribution develops into a log-normal 
distribution. This is further illustrated in the right 
part (e)-(h) showing the distribution of log(G) at the 
energy E = 0.01 eV just above the conduction band 
edge. 
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Figure 3: (a): Average resistance at E = 0.16 eV of uni- 
formly doped wires (circles) and surface doped wires (stars). 
The solid lines are estimates based on the single dopant 
transmissions (Fig. [TJ using J2). The MFP is found from 
the linear regions giving l e = 47 nm for the uniform distri- 
bution and l e = 190 nm for the surface doped wires, (b): 
ln(T) vs. L for the two different distributions in (a) (cir- 
cles and stars) and for uniform distribution at E = 0.01 eV 
(squares), (c): std(G) vs. L. Same symbols as above. 

The transmission, T = G ■ j^, decreases as T(L) ~ 
exp(— L/£) for wires in the localization regime. This 
is illustrated in Fig. [3] (b) showing (lnT) vs. L at 
E = 0.01 eV (squares) and at E = 0.16 eV for the 
uniformly doped wires (circles) and the surface doped 
wires (stars) . The localization length is obtained from 
the slopes in the final linear regions. At E = 0.01 eV 
we get £ = 7nm and for the uniformly doped wire 
at E = 0.16 eV we find £ = 133 nm. The sur- 
face doped wire does not reach a linear region in the 
considered length range and we can only say that 
£ > 200 nm. The calculated localization lengths gen- 
erally agree with the relation ([5]) giving £ = 8 and 
123 nm at E = 0.01 and 0.16 eV respectively. 

The length dependences of the sample-to-sample 
variations, std(G), shown in Fig. [3] (c) are qualita- 
tively different for the two radial dopant distributions. 
The purely surface doped wires (stars) approach the 
UCF level (horizontally dashed line) from below. On 
the other hand, the uniformly doped wires (circles) 
reaches a clear maximum before the UCF level is ap- 
proached from above. A similar maximum is observed 
at E = 0.01 eV (squares). Both maxima occur around 
L « ? e /2) which seems to be a quite general result 
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independent of energy and dopant concentration [9] . 
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Figure 4: MFP vs. energy for different radial dopant pro- 
files. The numbers in the legend refers to the fraction of 
the dopant atoms located at the surface. All solid lines are 
results obtained from the single-dopant transmissions. 



face doped wires as compared to more uniform dopant 
distributions. Moreover, the sample-to-sample varia- 
tions at wire lengths L < l e were smaller in the sur- 
face doped wires. These results can be understood and 
accurately estimated from the scattering properties of 
the single dopants, implying that relatively simple cal- 
culations are sufficient in practical nanoscale device 
modeling. 
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Using ((3|) we extract the MFP from the linear regions 
in Fig. [3] (a). By doing this at many energies, we 
obtain the energy dependent MFP shown in Fig. |4j 
The different markers correspond to different surface- 
dopant fractions. The solid lines arc obtained from the 
single-dopant calculations using (j4|) and correspond to 
the same surface dopant fraction as the closely lying 
markers. We see that the single-dopant estimates ac- 
curately resemble the sample-averaged results over the 
entire energy range. At all energies, the MFP increases 
with increasing surface dopant fraction, and at pure 
surface doping the MFP is on average more than three 
times larger than for uniform doping. The large in- 
crease in MFP for surface doped wires might be useful 
for device applications. However, there are potential 
problems with dopants close to the surface as they can 
be passivated by extra hydrogen atoms and thereby no 
longer provide extra carriers [5]. 

5. Conclusion 

By combining ideas of scaling theory and universal 
conductance fluctuations with DFT and Landauer for- 
malism we have analyzed the conductance properties 
of SiNWs. We have studied the transport in phos- 
phorus doped SiNWs by computing the sample aver- 
aged conductance, mean free path (MFP), localization 
length, and sample-to-sample variations as a function 
of Fermi energy, radial doping distribution, and wire 
length. The MFP was significantly increased in sur- 
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